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Abstract 

We investigate a holographic superconductor that admits an analytic treatment near the 
phase transition. In the dual 3+1 dimensional field theory, the phase transition occurs 
when a scalar operator of scaling dimension two gets a vacuum expectation value. We 
calculate current-current correlation functions along with the speed of second sound 
near the critical temperature. We also make some remarks about critical exponents. 
An analytic treatment is possible because an underlying Heun equation describing the 
zero mode of the phase transition has a polynomial solution. Amusingly the treatment 
here may generalize for an order parameter with any integer spin, and we propose a 
Lagrangian for a spin two holographic superconductor. 



1 Introduction 



Refs. [H [2] observed that much of the physics of superfluids and superconductors can be 
captured by a simple gravitational model. The connection, which relies on the AdS/CFT 
correspondence [3 [H [5], is that an instability for a charged black hole to develop scalar 
hair is dual to a superfluid or superconducting phase transition. Because the AdS/CFT 
correspondence maps classical gravity to a strongly interacting field theory, the observa- 
tion opens a window onto strongly interacting superconductors and superfluids where BCS 
theory and other weak coupling techniques are inadequate. 

The model of [H [2] has been studied over the last several years in great detail. See 
[8] for relevant reviews. The model, in its most basic incarnation, involves an abelian 
gauge field and a charged scalar $ in a space with a negative cosmological constant. 
Most of the studies have been numerical in nature. Solving the equations of motion that 
follow from the gravity action reduces to a set of nonlinear, coupled, ordinary differential 
equations which can then be solved on a computer. One notable exception is the p-wave 
superconductor studied in refs. [9j [10] where the dual 3+1 dimensional field theory has a 
vector rather than a scalar order parameter. The gravity model in this case involves a 
non-abelian gauge field and no charged scalar. The order parameter is a boundary value of 
an off-diagonal component of the gauge field itself. 

In all of these gravity constructions, the onset of the phase transition is triggered by a 
zero mode for the field that gives hair to the black hole. In the case of this 3+1 dimensional 
p-wave model, the zero mode can be solved for analytically, and with this solution in hand, 
one can study in detail properties of the field theory near the phase transition. 

It seems strange that an analytic treatment should exist for the more complicated p- 
wave model and not exist for the scalar model. In this paper, we show that an analytic 
treatment exists for the scalar model, in particular for a 3+1 dimensional field theory and 
a charged scalar whose mass saturate, the Breiteulonuer-FYeedmau bound [HlElE 

To a first approximation then, this paper is a rewriting of [10] for the charged scalar. 
We are able to compute the critical temperature T c of the phase transition. Then near the 
phase transition, we calculate the speed of second sound, the change in the free energy, 
and the growth of the order parameter. We can examine the phase diagram near T c as 
a function of superfluid velocity. We also calculate current-current correlation functions 
near the phase transition and study how the hydrodynamic poles move around the complex 
frequency plane. These Green's function computations offer independent confirmation of 
the sound speed calculated earlier in the paper. Finally, we discuss the conductivity and 
London equations for the field theory. 

We are also able to add a few embellishments in our rewriting of [lOj. Along lines 
suggested by jT5l [16] , we consider more general potential and kinetic energy terms for the 
charged scalar in the action, and we study how these modifications affect the phase transition 
and critical exponents. Without the modifications, the exponents are mean field. With these 
modifications, we find a class of holographic superconductors with critical exponents 7 = 1 
and 1 < 5 < 3. These results for the critical exponents do not appear to rely on having 

1 See however refs. |13l I14| for related analytic approaches to the scalar model. 
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an analytic expression for the zero mode or on the dimensionality of the systems The 
analytic expression for the zero mode allows us to find the constant of proportionality in 
the relations defining the critical exponents. In addition to finding the critical exponents, 
we are able then to see in detail how the second order phase transition may be spoiled by 
irrelevant terms in the gravity action. 

If an analytic holographic superconductor exists for scalar and vector order parameters, 
then perhaps it should exist for all integer spins. The differential equation for the zero mode 
in the vector and scalar case is a linear second order differential equation with four regular 
singular points, i.e. an example of a Heun equation [2D]. The special zero mode solutions 
we found turn out to be Heun polynomials. Moreover, the differential equation can indeed 
be generalized in a way that is suggestive of a tower of holographic superconductors with 
order parameters of integer spin. We conclude by writing down a Lagrangian for a spin 
two holographic superconductor coupled to a gauge field. There are well known difficulties 
involving causality, ghosts, and tachyons in constructing such Lagrangians, and we perform 
no consistency checks, leaving further investigation for future work [21]. An analytic spin 
two holographic superconductor is clearly an object worthy of study, hopefully bringing us 
closer to the real world of high T c superconductors with d-wave order parameters. 



2 The Gravity Dual of a Superfluid 

In this section, we consider a holographic model of a strongly interacting, relativistic, scale 
invariant superfluid [TJ [2j 122] [23] . The action for a Maxwell field and a charged complex 
scalar field coupled to gravity is 

S = 5bulk + Sbry (2.1) 

whera^ 



R + ^ | - -^F ab F ab -V-K\89- iAm 2 



2k 2 \ L 2 4e 2 " 



(2.2) 



and Sbry is a boundary action which ensures a well posed variational problem and also 
renders the on-shell action finite. We reveal Sbry below. For most of this paper, we will 
take 

A' = l + A|^| 2 , V = m 2 \^\ 2 + £|1'| 4 , (2.3) 

although in Section [2.41 we consider slightly more general formsQ 

This gravity action describes a dual strongly interacting field theory that undergoes a 
superfluid phase transition. If the bulk scalar field vanishes then the equations of motion 



2 While there exist earlier discussions of these critical exponents, our argument appears to be somewhat 
different. The previous literature includes refs. [171 118] where the authors make similar statements about 
the exponents in the mean field case. Refs. [151 1161 119] discuss how the critical exponents change with 
modifications to the action, but largely from a numerical perspective. 

''Roman indices a,b, . . . run from to 4 and are raised and lowered with the five dimensional metric g a b- 
The Greek indices n, v, . . . run from to 3 and are raised and lowered with the Minkowski tensor rf v with 
signature ( — h ++). 

4 We could also consider a nonminimal coupling of the form |^>| 2 |_F| 2 . This coupling shifts the mass of 
the scalar field in a way that spoils the analytic computation we will soon present, and so we tune it to zero 
by hand. 
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following from eq. ()2.2[) admit a black hole solution. Such a gravity solution corresponds to 
a thermal, symmetry unbroken phase of the field theory with a nonzero chemical potential. 
Following the general arguments in ref. pQ once the temperature of the black hole falls below 
a certain critical value T c , a new phase exists in which the scalar field condenses. A solution 
with a nontrivial profile for ^ in the bulk corresponds to a superfluid phase in the field 
theory. 

In this paper, we work in a weak gravity (or probe) limit n/eL <C 1 in which gravity 
decouples from the scalar and gauge fields (or Abelian-Higgs sector) [2J. In this limit, a 
solution to Einstein's equations is a black brane: 

ds 2 = ^ (-f(u)dt 2 + dx 2 + dy 2 + dz 2 + ^) , (2.4) 

where f(u) = 1 — (u/uh) 4 - At u = Uh > 0, there is a horizon, while in the limit u — > 0, 
the space asymptotically becomes anti-de Sitter. We will then study the phase transition 
involving the Abelian-Higgs sector in this fixed background space-time. The equations of 
motion in this limit are 

-±=(d-iA)aV=9K(B9-iA*) a = ^L + ^m-iA*\ 2 , (2.5) 



1 -_d c (V^F cd ) = e 2 g da J a , (2.6) 



y-9 

where the current is J a = iK[$>*(d - iA) a ^ - V(d + iA) a ^*]. 

From the usual rules of the AdS/CFT duality [H [5] the boundary value of the metric 
g^y acts as a source for the boundary theory stress tensor T^ v \ the bulk gauge field acts 
as a source for a conserved current corresponding to a global U(l) symmetry; and the 
near boundary data of the scalar ^> sources a scalar operator Oa with conformal scaling 
dimension A > 1. The conformal dimension of Oa is related to the five dimensional mass of 
the scalar field m 2 through m 2 L 2 = A(A — 4). In what follows, we will be considering the 
particular case of m 2 L 2 = — 4 which saturates the Breitenlohner-Freedman bound I12j. 
This mass yields an operator of dimension two. 

By boundary values, we mean more specifically and A^ in the following near 
boundary expansions of \& and A^: 

* = j^^ [h) u 2 Hu/5)-(0 2 )u 2 + ..) , (2.7) 

A, = 4 fe ) + ^<J> 2 + ... , (2.8) 

where we have introduced a UV cutoff 5 <C 1. Given these expansions, the one-point func- 
tions of the corresponding field theory operators can alternately be expressed as variations 
of the action with respect to the boundary values: 

(J") = lim - 1 , (CM = lim . 1 , (2.9) 
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where gfy 



we can 



lim u _ > o(^/ L) 2 g^u- Given the fact that (J^) and are canonically conjugate, 
identify with the chemical potential p and A^p with a superfluid velocity 
Having fixed the values of Ap and we need to add certain boundary counter- 

terms (Sbry to the action to have a well defined variational principle. There are the usual 
boundary cosmological constant and Gibbons-Hawking terms that are needed to regulate 
the Einstein-Hilbert action in an asymptotically AdS space-time. However, we take the 
probe limit and our focus is on the Abelian-Higgs sector so we ignore these terms. We do 
need counter-terms for the scalar: 



Sb 



ry 



d 4 x 



-7 (2\V\ 2 /L + (**n a d * + Vn a d a y*)) 



(2.10) 



Here 7 is the induced metric on a small radius slice u = 5 and n is a unit vector normal to 
the boundary and oriented outward. With these counter-terms one can check that in the 
limit -> the expansions (|2.7p and (|2.8p are consistent with the variational definitions 
(|2.9p although in the computation of higher point functions further counter-terms are often 
required. In the following, we find it convenient to introduce 



e = -V2(0 2 ) . 



(2.11) 



2.1 The Superfluid Phase 

In this section, we study the superfluid phase transition for a charged scalar with m 2 L 2 = —4 
in the probe limit. 

We begin by searching for solutions to the equations of motion that are independent of 
x, y, z and t. For convenience, we will set e = 1, L = 1 and Uh = ljfl We write the bulk 
scalar as \l/ = pe lip ■, and make a gauge transformation ^4^—7-^4^-1- db<p. In the new gauge, 
the phase 93 disappears from the equations of motion and the current becomes J a = Kp 2 A a . 
The n-component of Maxwell's equation (|2,6p now gives Kp 2 A u = which means we can 
choose A u = 0. The equations of motion reduce to 



K 1 / 2 u 3 d n 



u 6 



ud u [A' t /u] 
ud u [fA'Ju] 



A 2 

A 2 - -± 

1 f 



p(l + A P 2 ) + — — 
u z op 



KA t 



u 2 f 
^KAi , 



(2.12) 
(2.13) 
(2.14) 



where ' indicates a derivative with respect to u. At the horizon, the boundary conditions 
are that A t {uh) = 0, Ai(uh) < 00, p{uh) < 00. At u = 0, we fix the leading values of 



Here we define 



5 There is a minus sign discrepancy between this identification and the one in [22] 
Ai(0) = £j whereas in [22], £j was defined to be the phase gradient of the scalar. 

6 To restore Uh, quantities we introduce later such as the chemical potential fi, frequencies u, and wave- 
vectors k should be understood as the dimensionless combinations (MUh, ojuu, and ku^, respectively. To 
restore L and e, send ^/ — > tyeL, m — > mL, A — > A/e 2 L 2 , and B — > B/e 2 and multiply the Abelian-Higgs 
action by a factor of L/e 2 . 
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the vector potential -At(O) = p and Ai(0) = £j. We also set $W = 0. The normal phase 
corresponds to the solution = p(l — u 2 ), = 0, and p = 0. 

To search for the superfluid phase, let us begin by setting A4 = 0. We are interested in 
solutions where p is small: 

A t = cp + e 2 4> 2 + e 4 04 + 0(e 6 ) , (2.15) 
p = epi + e 3 p 3 + e 5 /05 + O(e 7 ) , (2.16) 

where e <C 10 The equation for At is solved at zeroth order by <fio = po (1 — u 2 ). In the 
special case po = 2 and m 2 = —4, the equation for p\ reduces to 



„ / /' 3\ , , ({l-u 2 ) 2 1 v 



and can be solved: 

Pi 



u 2 



1 + u 2 ' 

As described in Section HI this differential equation for pi is related to the Heun equation, 
an ordinary differential equation with four regular singular points. The existence of such 
a simple solution for p\ is guaranteed by the existence of Heun polynomial solutions for 
certain choices of the parameters in the Heun equation. 

Our strategy in solving for this power series in e will be to fix the fall-off e but to allow 
the chemical potential to be corrected order by order: p = 2 + e 2 5p 2 + e 4 5p,± + . . .. Thus, in 
solving the differential equations, we require the boundary condition that the 0(u 2 ) term 
in pi vanish (for i > 1) while 4>i(0) is allowed to be nonzero. The first few terms in the 
power series are 

2 = ^(l-n 2 )- ^ 1 "^ , (2.17) 



(l + n 2 ) 
n 4 

'24(l + u 2 ) 2 ' v ^ 48(1 + n 2 ) 



,3 = -(l-2A + B) nA / ^ +(1-8A-2B) U2 ^} 1+ ^ . (2.18) 



Note that 



^ = ^-(1 + 4^ + 45). (2.19) 
48 

If we invert the expression p = p c + 5p 2 e 2 + . . . , we find a result for the order parameter as 
a function of chemical potential near the phase transition: 



, V2 

5p 2 



(2.20) 



If we restore dimensions, then p should be replaced by puh = p/ttT. If 5p 2 > 0, then the 
new branch of solutions exists for p > p c or, fixing p and letting T vary, T < T c . Moreover, 
the order parameter grows as the square root of the reduced temperature, e ~ (T c — T) 1 / 2 . 
This critical exponent 1/2, often called (3, is a classic result from the Landau-Ginzburg 



7 A similar expansion was used in ref. [13] to extract an analytic expression for the coherence length and 
to derive the London equation near the phase transition. 
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mean field theory of phase transitions. We will be able to confirm below that the phase 
transition is second order by computing the free energy difference between the phases. 

The case 5p2 < is strange. The new branch of solutions exists only for T > T c , and 
we will see below, in the free energy computation, that the ordered phase has a higher 
free energy than the normal phase. Naively the ordered phase is unstable and the phase 
transition will not occur. We suspect that 5p2 < means one of two things. In the first case, 
the free energy difference will eventually become negative if we follow this new branch of 
solutions far enough, and T c corresponds to a spinodal point in a first order phase transition. 
In the second case, this new branch of solutions is a local maximum between the normal 
phase and a runaway direction in the potential landscape. At T c , the normal phase develops 
a runaway instability. The terms that cause 5p2 to become negative are indeed problematic. 
For B < 0, the p 4 term in the potential is inverted. Also, A < makes the kinetic term 
smaller. 

2.2 The free energy 

For the free energy calculation, we also need the $4 and p$ terms in the expansion of At and 
p. The expressions are too cumbersome to repeat here, but we do give the near boundary 
expansion: 

p = eu 2 + 0(u 4 ) , (2.21) 
A t = (2 + 5p 2 e 2 + 5p A e 4 + 0(e 6 )) - ( 2 + Q + <5/x 2 J e 2 + (2.22) 
-5 + 37 A + WB + 6(1 - 8A - 2B) log 2 



1152 



+ Sih e 4 + 0(e b ) u 2 + 0(u 4 ) , 



where 



5u 4 = - - ( - 1265 + 16,568,4 + 33205 - 62,96(L4 2 - 31.600AB - 4820S 2 
276,480 

+240(1 - 8A- 2B)(7 - 44A - 8B) log 2) . (2.23) 



The free energy of the field theory is determined by the value of the action (|2.ip evaluated 
on-shell, Q = —TS os . Employing the equations of motion, we may rewrite the bulk action 
(1231) (in the case A = B = 0) as 



d , x Vz^_ (rA(jA / + p(f) | u=a + _ ^ d uy/=gApA»f? 



(2.24) 



We have included a cut-off 5 to regulate the boundary divergences. For general A and B, 
the integral over u has more terms which we do not write. Using the explicit form of the 
metric fj2.4|) and the boundary counter-terms introduced above (|2.10p . the on-shell action 
reduces to 

1 



S os = 2p Vol 3 



pJ l + Zif + I H du ^ A^p 2 
Jo 



(2.25) 



With our homogenous ansatz for the fields, the integral over d 4 x yields a factor of V0I3/T 
where V0I3 is the spatial volume of the field theory. 
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For a background where = and 

A t = (2 + 5fi 2 e 2 + 5/M e 4 + 0(e 6 ))(l - u 2 ) , (2.26) 

the on-shell action is 

S vac = /3Vol 3 (4 + ^(l + 4A + 4 J B)e 2 (2.27) 

+ — - — (-1235 + 16,808,4 + 3560B - 62,480A 2 - 30,640 AB - 43405 2 
69,120 v 

+240(1 -8A- 2B){7 - UA - 8B) log 2)e 4 + <3(e 6 ) 

For the background where the scalar has condensed and hence ^ 0, we find in contrast 
that 

S* = /3Vol 3 (4 + ^(l + 4 J 4 + 4£)e 2 (2.28) 

+ 69 120 ^~ 1055 + 17 ' 258A + 42805 ~ 62,480A 2 - 30,640,45 - 4340B 2 
+240(1 - 8A- 25) (7 - UA - 8B) log 2)e 4 + 0(e 6 
Now 

AQ = n rf - fi va c = T(5 vac - S s{ ) = Vol 3 (-3^(1 + 4^ + 4S)e 4 + 0(e 6 )) . (2.29) 

Intriguingly, the sign of Ail, is correlated with the sign of 6/12 (|2.19p . Fixing /i and letting 
T vary, when the ordered phase exists for T > T c , the free energy difference is positive, and 
the ordered phase is less stable. When the ordered phase exists for T < T c , the free energy 
difference is negative, and the ordered phase is more stable. From the fact that the free 
energy difference scales as e 4 , we see that the phase transition is second order. For small e, 
e 4 ~ {p — /i c ) 2 - Thus, e 4 ~ (T c — T) 2 . The derivative of Aft with respect to temperature 
is continuous but non-differentiable. There exists a critical exponent a which describes the 
behavior of the specific heat C v = TdS/dT ~ (T c — T)~ a near a phase transition. From 
the behavior of the free energy, we see that a = 0, consistent with Landau-Ginzburg mean 
field theory. 

2.3 Superfluid Flow 

In this section, we generalize the background above to allow for the possibility of a superfluid 
flow. In terms of the bulk solution, this generalization requires turning on a constant value 
of Ai(0) = £j. Invoking rotational invariance, we take = (£, 0, 0) without loss of generality. 
We again solve the system in a small e expansion but we treat £ as a small parameter. In 
this double series expansion, we find that 

1 12 \ /1 „.2\ , 1 ,2 



( 2 + 2? J 0--*r) + <h<? + ... , (2.30) 

A x = e(l-^y £ 2 ) + ..., (2-31) 
" ^U..- (2-32) 



1 + u 2 s 4(1 + u 2 ) 
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These solutions can be used to compute the speed of second sound. In a two component 
fluid, there are typically two propagating collective modes, ordinary and second sound. In 
our probe approximation the two sounds decouple [22]; second sound is associated purely 
with the Abelian-Higgs sector while ordinary sound is associated only with the gravitational 
sector. 

The speed of second sound, like that of ordinary sound, can be computed from derivatives 
of the state variables. From ref. [22], the second sound speed squared in this probe limit 
should be 

2 _ dJ*/d£ 



2 dJt/dn 
The small u expansion for the vector field is 



(2.33) 



= 12-,.(7 + 4.4 + 4B) 

^ 1 + AA + 4B ' V ' 

4. = e- 6a /T 2 L ^ 2 + --- • (2-35) 

^ 1 + 4A + AB K ' 



Hence, up to higher order corrections in e, 

d~ « („-2)^ 1+ ^ + 4J? . (2.36) 

2 7 + AA + 4B ^ 1 8 7 + AA + 4B v ; 

This perturbative solution in £ can also be used to analyze the phase diagram of the 
system near the critical point /x c = 2. At the critical point, we expect the order parameter 
to vanish, so e = 0. The value of At at u = can be reinterpreted as the value of the 
chemical potential. These two facts give us a relation between the chemical potential and 
superfluid velocity along the critical line separating the two phases. We expect that 

+ ■ (2-37) 

2.4 More General Superconductors and Critical Exponents 

In this section, following a suggestion of [151 116j . we may consider a more general form for 
V and K. Consider 

K(p) = 1 + Ap a ~ 2 ; V(p) = l -m 2 p 2 + Bp a , (2.38) 

where 2 < a < 4. Then we need to consider a more general expansion for the fields of the 
form 

A t = p c (l ~ u 2 ) + e a ~ 2 a _ 2 + . . . , (2.39) 
</, = epi+e «-i pa _ 1 + ... . (2.40) 

We find that 4> a -2 = 5fi a -2{l — u 2 ), and hence 

. (2.41) 

OpLa-2 J 
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Thus, we can get any critical exponent /3 we desire between 1/2 and oo by tuning the value 
of a and making sure that 8p a ~2 > 0. As in the case 8p2 < 0, for bp a -2 < 0, we anticipate 
a first order phase transition or runaway pathologies. If we want a critical exponent less 
than 1/2, we first have to tune the values of quartic contributions such that Sp,2 = 0. Then 
we have the possibility of adding higher order terms with a > 4. 
In the particular case a = 3, we find 5p\ = (8A + 9B)/8 and 

p 2 = -(4^ + 3 B ) " 2 '°f; + 2 f . (2.42) 
8(1 + u z ) 

We have also calculated the free energy difference for this model, using the method detailed 
in section 12.21 The result is 



A Q = ft sf - ft vac = Vol 3 (--^ (8,4(5 -6 log 2) + 95(3-4 log 2)) e 3 + 0(e 4 )) . (2.43) 

We see that Afl ~ (T c — T) 3 and thus the critical exponent a = —1. 

So far, we have focused on two critical exponents, a and j3, describing the behavior 
of the specific heat and the order parameter near the phase transition. There are several 
more, although general scaling arguments demonstrate that only two of these exponents 
are independent. The rest are fixed through scaling relations such as a = 2 - /3(J + 1) and 
7 = f3(5 — 1). The exponent 5 is defined at T = T c through the relation 

$r( 6 ) „ (0 2 ) 5 (2.44) 

while 7 is the susceptibility for T <T C : 

d(0 2 



(T c - Tji . (2.45) 



As we have seen above, f3 is fixed by the form of the series expansion (|2.39|) and (|2.4U[) . 
and in particular by the subleading term in At. The critical exponent 5 is fixed in a similar 
fashion, by the subleading term in p. In the previous examples, we solved for p a -i by 
insisting on three conditions: that p a -i be finite at the black hole horizon and that the 
leading and subleading terms at the boundary vanish. In general, we would not be able 
to impose three conditions on a second order linear differential equation, but we were able 
to tune the value of 5p a -2 that we had left arbitrary in solving for (j) a -2- We could have 
solved this system in a slightly different way, setting <5/i a -2 = 0, and thus freezing p,/T to 
the critical value. Still imposing regularity at the horizon, we would find that we need to 
turn on a source $W for the operator O2 in the solution for p a -\- Hence, there would be a 
scaling relation of the form (|2.44|) where 5 = a — 1. 

For the class of models we are considering, we see that (5 = l/(a — 2) and 5 = a — 1. 
From the scaling relations, the susceptibility will always be 7 = 1 while the heat capacity 
is described by 

o = — ■ (2-46) 
a — I 

This value for a agrees with our two free energy calculations in the cases a = 3 and a = 4. 
It would be interesting to verify this relation explicitly for all a such that 2 < a < 4, but 
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we have no reason to doubt its veracity. These results appear to be independent of the 
space-time dimension and of whether we have an explicit analytic solution for p\. What 
knowledge of p\ provides are the constants of proportionality in the definitions of the critical 
exponents and confidence that the series expansion in e is well behaved. 

Previous papers discussing critical exponents in a similar context are [151 EH EH OS CES] • 

3 Green's Functions 

To calculate the Fourier transformed retarded current-current correlation functions, we need 
to study fluctuations of the gauge and scalar field in our background. Our background is 
the black D3-brane (j2.4j) with chemical potential fi > ji c close to the critical value and no 
superfluid velocity. (There is no obstruction to introducing a small nonzero £ 7^ but the 
analysis becomes more complicated.) 

We will consider two types of fluctuations, a transverse fluctuation 5A X ~ e —i&t+iky 
and a longitudinal fluctuation 5A X ~ e - l ^ t + tkx _ The first type of fluctuation contains no 
propagating modes while the second will allow us to look at the second sound and shear 
modes. To keep the analysis simple, we restrict to the hydrodynamic regime where the 
order parameter, frequency, and wave-vector are small compared to the temperature. In 
our dimensionless notation, e, k, oj <S 1. We will also set A = B = in the following, 
choosing the simplest possible forms for V and K ()2.3|) . 

3.1 Pure Transverse Mode 

The transverse mode is described by fluctuations of the field A x with only y spatial depen- 
dence. We decompose the fluctuations into Fourier modes: 

5A x (u,t,y) = a x (u)e^ t+ik y . (3.47) 

These fluctuations decouple from the other scalar and gauge field components and satisfy 
the second order differential equation, 

a ' + {j-u) a ' + {— — 57 

Near the horizon u = 1, we find that a x ~ (1— u)^ llu '^ satisfies either ingoing or outgoing 
plane wave type boundary conditions. Consistent with the presence of an event horizon, it 
is natural to choose ingoing boundary conditions (the minus sign in the exponent). This 
choice leads to retarded, as opposed to advanced, Green's functions in the dual field theory 
[24] . At the boundary u = of AdS, we would like the freedom to set a x (0) = a x o to some 
arbitrary value of our choosing, corresponding to perturbing the dual field theory by a small 
external field strength. These two boundary conditions along with the differential equation 
uniquely specify the functional form of a x . 

While an analytic solution of (j3.48|) does not appear to be available, one can easily solve 
this equation in the limit of small u, k, and e. The first few terms in this expansion are 

/ 1 _ u 2 \ 

a x = a x0 I — — — 2 ) (l + <?a xt + e 2 cja X0Je + k 2 a xk + io 2 a XU) + . . .) . (3.49) 



(3.48) 
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We find 

,2 



WT^) • <3 ' 50) 



iu 2 



(3.51) 



16(1 + u 2 ) ' 

a xk = i ln(u) In - ^Li 2 (u 2 ) + ^Li 2 (u 4 ) . (3.52) 

The expression for a XUJ is too cumbersome to give here. Near the boundary, this solution 
(|3.49p has the expansion, 

fiu e 2 we 2 u? log 2 1 2 2 / 1 \\ , 

Ox = axo + Oxo I y - -g - Yg 4 + 2^ w ~ >\2~ ^ J Ju +... . (3.53) 

From this near boundary expansion, we can calculate the two-point function for the current 
in the hydrodynamic limit: 

, , _ I 'iu e 2 iuje 2 w 2 log2 9 l9 , \ , nrl , 
GP B (w,* ! ) = 2fy---— ^ + ( W 2 -/c 2 ) C J . (3.54) 

This Green's function is morally a second derivative of the on-shell action with respect to 
a X Q. Note that there is a counter-term ambiguity, proportional to an arbitrary constant c, 
of the form discussed in [TOj [25] . 

3.2 Second Sound Mode 

In general, second sound modes are expected to produce poles in the density-density cor- 
relation function. We thus need to consider fluctuations in the conjugate field At. A 
self-consistent set of fluctuation equations also involves the scalar field and longitudinal 
fluctuations of A x : 

SA t (u,t,x) = a t {u)e~ iult+ikx , (3.55) 
5A x (u,t,x) = a x (u)e- iut+ikx , (3.56) 
6V(u,t,x) = if;(u)e- iult+ikx /V2 . (3.57) 

These four fluctuating modes satisfy four second order differential equations (if one also 
counts the equation for ip*) and one first order constraint equation: 

a'l--a' t - + %-^a x - -^- {^ + 2^-^-2^*) = 0,(3.58) 
u V u J J J 2u J 

(f 1\ /a; 2 p 2 \ ku kpty-r) n , q , Q , 

ax+ {j-u) ax+ {p-^f) ax+ p at+ 2u 2 f = °' (3 - 59) 



/' 3\ [ (co + cj)) 2 A-k 2 u 2 \ (a, + 2<ft)p kp 

l~u)^ + { p + u 2 f )^ + — J 2 — at + j a - = Ms.eo) 

ua' t + kfa' x + ((^* - ij))pl - (if;*' - i//)p) = 0.(3.61) 
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We have checked that the derivative of the constraint equation (|3.6ip with respect to u 
is a linear combination of all five differential equations (|3.58p - (|3.6ip . Thus if a solution 
of the first four differential equations satisfies the constraint for some u, it will satisfy the 
constraint for all u. 

There are seven integration constants associated with this linear system f|3.58j) — (|3.61j) . 
If we look at the horizon of the black hole at u = 1, we find seven different kinds of behavior. 
There exist six solutions that have plane wave behavior for a x , ip, and ip* near the horizon 
of the form (1 — u) 1 * 1 "^/ 4 . There is also a pure gauge solution, 

at = —iui , a x = ik , ip = ip , ip* = —ip . (3.62) 

As in the analysis for the transverse fluctuations, we choose pure ingoing boundary condi- 
tions corresponding to (1 — u)~ lul / 4: behavior. At the boundary u = of our asymptotically 
AdS space, we would like to be able to perturb the system with arbitrary boundary values 
of at and a x but set the source terms for ip and ip* to zero. These are four constraints and 
we have only three ingoing solutions. Thus we will also need to make use of the pure gauge 
solution to enforce the boundary conditions at u = 0. 

We solved this system perturbatively in oj, k, and e. We present the results here in the 
limit where oj ~ k 2 ~ e 2 . The near boundary expansion (u = 0) of the solution takes the 
form 

at0 + k ^ + ax0 u ^ . fc4 + ? . g4 + M8e2w _ 48Q . w2 

+8k 2 (5ie 2 + 12w)] u 2 + ... , (3.63) 
axQ _ ^tok + a^u: ^ . fc4 + ^ + _ m ^ 2 

+8k 2 (5ie 2 + 12w)] u 2 + . . . , (3.64) 



a t 



ip = _^8k[l2/c 4 + (7e 2 -120io;)w + 3A: 2 (e 2 + 16(l-i)a;)]u 
-^-iek [48/k 4 + e 4 + 12(2 - 3i)e 2 oj - 96(3 + 4»> 



^8ie [12A; 4 + (7e 2 

a *0^urAOlA , £ 4 + 12(2 _ 3 . )e 2 w _ 96(3 + 4 - )w 2 

+16A: 2 (e 2 + 6(1 - 2i)u))] u 2 + . . . , (3.65) 

where the pole structure is given by 

V = 960w 3 + 56i(12A; 2 + e 2 )uj 2 - 12(16A; 2 + 3e 2 )k 2 oj 

-i(48k 4 + 16A; 2 e 2 + e 4 )A: 2 . (3.66) 

Note the expression atok + a X QCO is not homogenous in our scaling limit. We leave this piece 
in the Green's functions to make it obvious that the Green's functions satisfy the Ward 
identity k fJ/ G^ u = 0. We have included the leading corrections proportional to ato and a X Q. 
There are terms in the expansion proportional to u 2 In u that give rise to the counter-term 
subtraction ambiguity mentioned above but they are subleading in oj, k, and e. 

Let us study the pole structure in two different limits. In the limit k <C e, the poles are 
located at 

lie 2 89ik 2 

" = "120 "W < 3 ' 68) 
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The first two poles are propagating modes that we identify with second sound. Indeed, the 
speed of second sound agrees with the earlier result (|2.36[) from Section \2. 31 The position of 
the third pole is determined mostly by the size of the order parameter e and so we associate 
it with the zero mode that causes the phase transition from the superfluid phase back to 
the normal phase. 

In the opposite limit, k 3> e, where the order parameter is small, the behavior should 
be close to that of the normal fluid. In this limit, we find 

±3 — z 2 ±9-i 2 

U = ~^~ k + ^40- £ ' (3 - 69) 

ik 2 ie 2 , . 

- = — "20- < 3J0 > 

The first two poles are associated with the zero modes that cause the phase transition from 
the normal phase to the superfluid phase, while the third pole is associated to a diffusive 
mode of the conserved charge density. Indeed the location of this diffusive pole is determined 
by the dynamics of the normal phase and was calculated without the order e 2 correction, 
long ago [26]. As we vary e and k, the number of poles cannot change. The two zero mode 
poles evolve into the sound poles of the previous limit while the diffusive pole becomes the 
zero mode pole of the previous limit. 

As a further check, we consider a particular static limit of the density-density correlation 
function. From (j3.63|) . we can read off the Green's function, 

k 2 

G tt (u), k) = -2— (A8ik 4 + 7ie 4 + 148e 2 w - 480iw 2 + 8k 2 (5ie 2 + 12a;)) . (3.71) 

We are interested in the long wave-length limit of this Green's function: 

UmG tt (0,k) = 14 . (3.72) 

This long wave-length limit is equal to a thermodynamic susceptibility, 

r) 2 P r)P 

HmG tt (0,fc) = ^ = ^- . (3.73) 
Given this relation, we see that (I2.34p agrees with (I3.72p . 

3.3 Ohm's Law and London Equations 

Via Ohm's Law, a homogenous, time dependent electric field 5A X ~ e - tult should produce 
a current in the system. To investigate this long wave-length limit of the conductivity, we 
set k = for the two-point functions computed above. We have 

r( W ,0) = 2^ +! ^-^ W + CW 2 + O( W 3 )j . (3.74) 

Note that the k — > limits of G xx computed in Sections 13.11 and 13.21 agree. 

Identifying the electric field E x = iuj5A x and recalling Ohm's Law, the conductivity is 
related via a Kubo type formula to the retarded Green's function 

<T*x(u>) = ■ (3.75) 

ILU 



13 



The term proportional to e 2 in G xx thus produces a pole in the imaginary part of the conduc- 
tivity. As discussed in [21 [27] , by the Kramers-Kronig relations (or by properly regularizing 
the pole), there must be a delta function in the real part of the conductivity, indicating 
the material loses all resistance to DC currents and suggesting the phase transition is to a 
superconducting state. Here we can calculate the residue of the pole analytically close to 
the phase transition: 

ie 2 

Res w=0 a xx = 2— . (3.76) 



An important observation is that in our system, the oj — > and k — > limits of the 
Green's functions commute. The residue of the pole in the conductivity is related to the 
long wavelength limit of the current-current correlation function in the following way: 

iRes^ocJxx = lim lim G xx (lo, k) . (3.77) 
The limit in the opposite order is related to a thermodynamic susceptibility: 

where £ is the superfluid velocity. It follows from (|2.35p that 



8 2 P 8T X e 2 

w=~k=- 2 j- < 3 - 79) 

As emphasized in this context in [27] , that the limits commute implies the system really 
does become a superconductor below T c . Given that the limits commute, the system obeys 
a London type equation for k and uj: 

(J x ) = -2^A X . (3.80) 

o 

If we now imagine that the U(l) symmetry is weakly gauged, these London equations imply 
not only infinite DC conductivity but also a Meissner effect with a London penetration 
depth that scales as A ~ 1/e ~ (T c — T) 1 / 2 . 

4 Heun Equation 

Consider a charged spin-s field h aia2 ... as propagating in a curved space-time of the form 
(|2.4p in d + 1 dimensions, 

ds 2 = ^ ( -f(u)dt 2 + dx 2 + ^f- ) , (4.81) 



% (-f(u)dt 2 + dx 2 + 
u 2 V /(«) 

where x = (xi,X2, • • • ,x^_i). We turn on a nonzero radial electric field generated by the 
potential At(u), and an off-diagonal, purely spatial component of h of the form fy(u) = 
hi 1 i 2 ...i s where i indexes x only. Putting aside the well known difficulties of constructing 
actions for charged, higher spin fields in curved space-time, let us suppose for the moment 
that this component ^ satisfies an equation of motion of the form 

1 



-9 



d u {V^gg ilh g i2i2 ■ • ■ g isis g uu d u t?) - g hil g i2i2 ■ ■ ■ g isia {g tt A 2 ^ + m 2 ^) = . (4.82) 
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Given our ansatz, this equation of motion can be written in the form 

*» + + *» + (f - ^ ) *(«) = • (4-83) 

This equation governs the onset of a superconducting instability in two different known 
holographic systems. In both cases, the holographic system is in the so-called probe limit 
where gravity is weak and the background geometry is fixed. In the case s = 1 and m = 0, 
the equation governs the instability in the non-abelian superconductor of [281 129] 130] . In the 
case s = 0, the equation governs the phase transition in the s-wave superconductor of [HE]. 
We will see in the next subsection that (|4.83p may arise from a charged spin-2 field as well. 
In the original papers [TJ [2j [28], the focus was on field theories in d = 2 + 1 dimensions, 
perhaps because of experimental relevance for high temperature superconductors which are 
typically layered materials. We will choose .At (it) = fi (1 — u 2 ) and f(u) = 1 — -u 4 , implying 
the field theory is d = 3 + 1 dimensional. 

In d = 3 + 1 dimensions, the authors of [9] noted that ([4.83P has an analytic solution 
for [i = 4, s = 1, and m = 0. This observation was later developed by jlOj to calculate 
transport properties of the nonabelian superconductor near the phase transition. In fact 
[10] noted the existence of a tower of polynomial solutions for /i = 4n, with n a positive 
integer. 

Closer scrutiny of ([4.83p reveals that these special solutions are actually Heun polyno- 
mials [20J. The Heun differential equation has the form 

d 2 y ( 7 5 e \ dy af3x 



dx 2 \x x — 1 x — a) dx x(x — l)(x — a)^ ' ^ ^ 

where a + /3 + 1 = 7 + S + e. The differential equation ([4.83P can be transformed into ([4.84P 
with a = — 1 with the identifications 

,2 



m 2 L 2 = 4(s — 1), // = 2(/c -(- s + 1) and x = u 2 . One finds that 

a = —k , /3 = —k — s , 7 = s + 1 , (5 = 1, e = — 1 — 2A; — 2s , 

and g = A;(fe + s+ 1). There exist polynomial solutions of this equation for all non-negative 
even integers k. The first few polynomials are 

2/0 = 1, (4.85) 

y2 = l- 2 Jl±^ x + x \ (4.86) 
1 + s 

4(5 + s) 2(54 + 25s + 3s 2 ) 2 4(5 + s) , 4 

2/4 = 1- '-x + —^x 2 - — -x 3 + x 4 . 4.87 

y 1 + 8 (1 + S)(2 + S) 1 + 8 V ' 

4.1 Spin-2 fields 

The Pauli-Fierz Lagrangian for a spin two field in flat space takes the form 

£pf = \ [h ab > c h ab , c - 2h a h a + 2h a h a - h a h a - m 2 (h ab h ab - /i 2 )] , (4.88) 
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where h a = h ab b and h = h a a . In making this action generally covariant, ambiguities arise 
because of possible couplings to the Riemann curvature tensor: R abc dh ac h bd , R ab h ac h b c , 
R ab h ab h, Rh ab h ab , and Rh? . Note that while h c h c is equivalent to h aC:b h bc,a in flat space, in 
curved space, a lengthy integration by parts demonstrates that 

(V a h ac )(V b h b c ) ~ (V b h ac )(V a h b c ) + R ab h ac h b c - R c dba h ad h b c . (4.89) 

For our spin two field coupled to a U(l) photon, we take our action to be of the form 

(D c h ab Y(D c h ab ) - 2(D a h ac )*(D b h b c ) + [(D b h ba f{D a h) + c.c] - (D a h)*(D a h) 
+2R abcd {h a yh bd + 2R ab {h a yh b c -m\{h ab Th ab -h*h^ , (4.90) 



where D = V — iA. While this action is manifestly gauge invariant, it is not clear whether 
it is free of ghosts and tachyons. 

We fix the only nonzero components of h ab to be h xy = h yx . We take the background 
to have the form (|2.4p and we assume that the only nonzero component of the vector field 
is A t . Assuming a homogenous ansatz where h xy depends only on the radial coordinate, u, 
the resulting equation of motion is 



h" J r[- + -T)h'+[-±- -^y- )h xy = 0. (4.91) 



1 f\ , (A\ _ m*L 2 

With the choice m 2 L 2 = 4 and At = fi(l — u 2 ), this equation will have Heun polynomial 
solutions for discrete values of [l of the form discussed above for s = 2|f| 

5 Discussion 

The Abelian-Higgs model in d = 3 + 1 dimensions with m 2 L 2 = —4 we considered here was 
studied before numerically in ref. [32] as part of a survey of different models of holographic 
superconductivity. We can confirm several of the numerical results in their Table 2. The 
confirmation requires a little work in order to restore units and to present the results in the 
same statistical ensemble. In the canonical ensemble (fixed charge density), we find from 
He = 2ttT, (HJZOD, and (IBTTil) that 

{Oi, » -{^-* 2 T!I\T C -Tf\ (5.93) 

72 

limG xx (w,0) w -—tt 2 T c (T c -T). (5.94) 
ui— >o 7 

We have left e = L = 1 in the above expressions and set A = B = 0. Up to a factor of 
two difference in the definition in J* and a minus one in the definition of (O2), these results 
agree well with Table 2 of [32]. 



While this work was being prepared, ref. [31] appeared with a discussion of d-wave holographic super- 
conductors for a different choice of action. 
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One interesting aspect of our results is their dependence on the modifications to the 
kinetic and potential energy of the charged scalar, parametrized by the constants A and B 
and the exponent a. The growth of the order parameter, the speed of second sound, and the 
behavior of the free energy near the phase transition are all sensitive to A and B. We saw 
that for K and V of the form (|2.3p . if 1 + AA + 4B < 0, the second order phase transition 
is destroyed. Given the existence of a second order phase transition, the critical exponents 
do not depend directly on the values of A and B, but they were sensitive to the choice of a 
in (|2.38p . While modifications to the action of the form (|2.3|) seem to be generic in stringy 
embeddings in the large N and large coupling limit, it seems likely that modifications of 
the form (|2.38p will be suppressed by inverse powers of N and the coupling strength^ 

The sensitivity of our results near T c on A, B, and a suggests that as we go to lower tem- 
perature, the properties of the system will depend more and more on higher order, irrelevant 
terms in the gravity action. This dependence makes our approach to superconductivity both 
richer and riskier. Our task is richer because at the phenomenological level there clearly 
exists a wealth of gravity theories to be explored and mined. Our task is riskier because it 
is less clear that all of these gravity theories make physical sense interpreted as dual field 
theories. The most well established AdS/CFT dualities are between ten dimensional string 
theories and supersymmetric field theories living in 3+1 space-time dimensions. Without a 
stringy embedding of (|2.2p . it is difficult to describe exactly the dual field theory or even, in 
fact, to be sure that such a dual exists. While a generic gravity dual will possess quadratic 
terms in the action of the form ()2.2|) . higher order terms will depend on the details of the 
stringy embedding. 
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